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Abstract 

We discuss tests of the charged Higgs sector of the Two Higgs Doublet Model (2HDM) of Type 
II in the light of recent flavour physics data. Particular attention is paid to recent measurement of 
purely leptonic decays of heavy-light mesons, which depart more or less significantly from the Standard 
Model (SM) predictions. We derive constraints on the parameters of the 2HDM type II from leptonic 
and semileptonic AF = 1 decays as well as loop processes (b — > sj, BB mixing or Z — > bb) sensitive 
to charged Higgs contributions. 

The outcome of this work is that while 2HDM Type II can fit individual observable through 
fine-tuning schemes, in a combined analysis it does not perform better than the SM by favouring a 
decoupling solution. Assuming that 2HDM Type II is realized in Nature, constraints on its parameters 
(m H + and tan/3) are derived. A limit on the charged Higgs mass m H + > 316 GeV at 95% CL is 
obtained irrespective of the value of tan (3. This limit is dominated by the 6—7-57 branching ratio 
measurement. 

All results have been obtained with the CKMfitter analysis package, featuring the frequentist 
statistical approach iffit to handle theoretical uncertainties. 
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1 Introduction 



The overall agreement between the Standard Model (SM) and data in the quark sector is 
particularly impressive: Flavour- Changing Neutral Currents (FCNC) are as small as predicted 
in the SM and the Kobayashi-Maskawa (KM) mechanism has been proven to describe the 
observed CP-violating phenomena in flavour physics with a good accuracy. One example of 
this convergence is provided by the global fit of the CKM matrix elements [I] illustrated in 
Figure [TJ All constraints (either loop and tree observables or CP- violating and CP-conserving 
quantities) point towards a unique solution, which proves the CKM mechanism to be at work 
in flavour transitions (within the present accuracy) and establishes the KM mechanism as a 
dominant source of CP violation in K- and 5-meson systems. The consistency between these 
predictions is deeply related to the fact that a single Higgs doublet provides a mass to all 
fermions through Yukawa couplings in the SM and that all but one CP-violating phases can be 
rotated away by redefinition of the fields. 
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Figure 1: Superimposed individual constraints at 95% CL for the SM global fit. The yellow bean is 
the solution driven by the combination of all individual constraints at 95% CL. Regions outside the 
red circle are excluded at 95% CL. 

Extensions of the SM are often based on the introduction of additional fields interacting with 
the quarks (gauge bosons of new interactions, supersymmetric particles, technifermions [2|. . . ). 
These new fields bring along new arbitrary parameters, often inducing dangerous FCNC pro- 
cesses, as well as new CP-violating phases that cannot generally be rotated away. Therefore, 
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flavour physics and CP-violation data provide stringent constraints on the parameters of these 
extensions, which yield generally fine-tuning problems. From this point a view, a rather min- 
imal extension of the SM, with a limited number of new parameters to fix, consists in the 
two-Higgs doublets models (2HDM). Indeed, in the SM, the same doublet is used to couple the 
left-handed fermion doublets with right-handed up- and down-type quarks at the same time 
(exploiting the fact that the representation of the Higgs doublet under SU(2) is pseudoreal, 
and that the hypercharges of the scalar doublets coupled to right-handed up- and down-type 
quarks are opposite). This choice is imposed by no other arguments but economy, and one 
can introduce two complex doublets 0i and 02 of opposite hypercharge, rather than a single 
doublet [21 IU El El E] • One assumes that electroweak symmetry breaking occurs because 
the neutral components of the two doublets acquire two a priori different vacuum expectation 
values: 

^ l/nU_ln\l _ ( V ^/V2 \ 

,2 i „,2 _ „,2 



KOK|o>[=( ..J^), M<h\o)\ = { j 'q ), (i) 



denoted respectively V\ and v 2 (with the constraint v\ + v\ = v , where v is the SM Higgs 
vacuum expectation value). 

2HDM models contain 8 degrees of freedom in the Higgs sector, out of which 3 are used to 
provide a longitudinal polarisation to the weak gauge bosons W and Z. Five real fields remain: 
two charged Higgs fields H , a neutral pseudoscalar Higgs A and two neutral scalar Higgs (h° 
and H°). The additional parameters needed to describe this SM extension are the masses of 
H ± , H° and A, the ratio of vacuum expectation values tan/3 = V2/V1 and an angle describing 
the mixing between h° and H°. 

Different versions of the 2HDM were labeled according to the couplings of the Higgs doublets 
to the quarks [6|: type I corresponds to 0i coupling to both up- and down- type quarks whereas 
02 does not couple to any quark, type II corresponds to 0i coupling to down-type quarks 
whereas 02 couples to up-type quarks (and to leptons), and type III to 0i and 02 both coupling 
to both types of quarks. Among these various possibilities 2HDM Type II is particularly 
alluring, because of its resemblance with the SM in the quark sector. One has two Yukawa 
matrices y d,u describing the couplings among quarks (and one for the lepton sector, y e ): 

C n ,Y = -Q L <Piy d D R - Q L <p 2 y u U R - L L <p 2 y e E R + h.c, (2) 

where Ql and Ll denote left-handed fermion doublets, and D R , U R , E r down-, up-type and 
charged-lepton right-handed singlets, defined in a similar way as in the SM (actually, the SM 
is recovered through the identification 02 = iu 2 (t>\- l with a 2 the complex Pauli matrix). 

One has to re-express these couplings in terms of mass eigenstates. The structure of the 
Yukawa terms yields a SM-like structure for the quark sector: there is a CKM matrix which 
is the only source of flavour-changing interactions and there are no flavour-changing neutral 
currents at tree level. But there are new flavour- changing charged interactions, corresponding 
to the exchange of a charged Higgs rather than a W (obviously, there are also interactions of 
quarks with neutral Higgs fields, as well as couplings of the Higgs fields to the leptons). Once 
quark and Higgs fields are expressed in terms of mass eigenstates, one obtains the following 
charged-Higgs interactions for quarks and leptons: 



" r 

ij 



tan /3^-ULiVijdnj + cot P^-u Rj V i:j d Lj + tan/3^-u Lj £ Rj 
M w M w M w 



H + + h.c. 

(3) 
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Therefore, the 2HDM Type II provides a very interesting extension of the SM: it exhibits 
the same CKM structure, has a natural mechanism of suppression of FCNC plaguing many 
other models, but it exhibits a different structure for charged currents by the addition of new 
(scalar and pseudoscalar) interactions. Furthermore, it relies on a limited number of additional 
parameters, i.e., the mass m H + of the charged Higgs and the ratio tan/3 of the couplings to 
up-like over down-like quarks (if we restrict our study to charged currents). In many situations, 
the change induced by the 2HDM Type II amounts to a redefinition of some of the parameters 
already occurring in the SM expressions, with a new dependence on and tan/3. 

Eventually, in addition to that predictability virtue, the 2HDM Type II is embedded into 
the most simple supersymmetric extensions of the SM (MSSM), at least at tree level jl] (for 
large tan (3, loop effects might lead supersymmetric theories to coincide with a 2HDM Type 
III rather than the 2HDM Type II). Searches for such charged Higgs are obviously among the 
prospects of the LHC experiments ATLAS and CMS [§]. 

Since the decays mediated by a weak charged current are an excellent laboratory to search 
for charged Higgs boson contributing in addition to W ± bosons, we have collected the measured 
decays potentially sensitive to contributions from charged Higgs, for which a good control of 
the theoretical hadronic uncertainties can be achieved. A combined analysis of their branching 
ratios in the light of the 2HDM Type II is then performed within the frequentist statistical 
scheme developed by the CKMfitter group [TU] . 

It is convenient to categorize these observables as follows: 

1. The leptonic decays of mesons mediated by quark-annihilation at tree level T[K — > 
Iiu]/Y[k — > fj,u], B[D — > fiu], B[D S — > /xu], B[D S — > ru] and B[B — > tu], where B stands 
for branching ratio and T for the decay width. In addition we also consider the strange 
hadronic decay of the r lepton, r — > Ku, through the ratio T[t — > Ku]/T[r — > tiu], which 
can be seen as reversed leptonic decays. 

2. The semileptonic decays B — > Dtu, through the ratio B[B — > Dtu]/B[B — > D£v], and 
K — > irlv, through the ratio B[K — » 7tLiu]/B[K — > weu]. 

3. The B® and B® oscillation frequencies, Am^ and Am s . 

4. The Z partial width into b quarks, Rf, = Y[Z — > bb]/T[Z — > hadrons], which exhibits 
electroweak charged currents through Zbb vertex radiative corrections. 

5. The measurement of the FCNC radiative decay b — > 57 through the ratio B[B — > 
X sl ]/B[B^X c ev\. 

Most of these observables, either from tree (first two categories) or loop (last three ones) 
contributions, are established individual benchmarks to constrain or measure the parameter 
space of the 2HDM Type II or cognate supersymmetric models [TTl [T2"] . The measurements 
of these observables and their experimental uncertainties are displayed in Figure [2] together 
with their SM prediction at 95% CL. Figure [2] shows an overall fair agreement between the 
various observables and their SM predictions, with the notable exception of the tauonic B + 
decay B — > tv (the deviation is yet lower than 3 cr). Let us notice that the recent CLEO 
measurements [13j of D s —t ilv and D s — > tv decays are now in agreement with their SM 
predictions in contrast to the situation reported at the time of the 2008 summer conferences |H3 . 

One of the main objectives of the analysis reported in this article consists in determining 
whether a 2HDM Type II can accommodate the discrepancy coming from the large value of 
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Figure 2: Comparison of the measurements relevant to constrain the 2HDM Type II and their predicted 
value in the SM. The black dots indicate the deviation of the experimental result from its prediction, 
assuming Gaussian distributed errors for both theory and experiment. The deviation is expressed as a 
signed significance whith positive values indicating a mesurement higher than its prediction. Note that 
this figure is for illustration purposes only, since all the errors are treated as Gaussian here. In the 
rest of the paper, we use the ~Rfit prescription to deal separately with statistical and systematic errors 



B —7- tv as measured by the B factories. A second aim is to determine the allowed region 
of the parameter space (m H , tan/3) as constrained by above set of low-energy observables, 
corresponding to AF = 1 tree processes or loop-induced processes featuring a single charged 
Higgs exchange. This region of parameter space can be compared to the limits set by LEP 
from the (absence of) direct production of charged Higgs bosons. 

2 Observables and theoretical context 

This section details the theoretical predictions for our set of observables in the SM and how they 
are modified in the context of charged Higgs contributions. A summary of the relevant mea- 
surements and parameters together with their corresponding uncertainties is given in Tables [1] 
andd 

2.1 Leptonic decays 

The decay of a charged meson M into a leptonic pair £u£ is mediated in the SM by a charged 
weak boson, with the branching ratio: 

B[M -+ Ivthu = (l - ^f) 2 \V quqd \ 2 f M r M (l + 6»£), (4) 

where q u (qa) stands for the up-like (down-like) valence quark of the meson respectively, V quqd 
is the relevant CKM matrix element, /a/ is the decay constant of the meson M (describing 
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Input 



Value Unit Accuracy Reference 



Branching Ratios 



T[K -»- fw]/T[iT \iv\ 


1.336 ±0.003 






(0.2%) 


[15j 


T[K° -»■ 7r/iz/]/r[^° -»■ ei/] 


0.6640 ± 0.0026 






(0.4%) 


m\ 


r[r -»■ ^]/r[r -»■ 7ri/] 


6.370 ±0.215 


10" 


-2 


(3.4%) 


[16] 


->■ H 


3.82 ±0.32 ±0.09 


10" 


-4 


(8.4%, 2.4%) 


m 


io r 7 ^» , i 

d[.D s — > /Ltl/J 


r Aft | r\ A r\ 

5.93 ± 0.40 


10" 


_3 


(6.7%) 


[18] [13] 


£[D S -> rz/] 


5.62 ± 0.44 


10" 


-2 


(7.8%) 


Lisj Li3j 


B[S -> ri/] 


1.73 ±0.35 


io- 


-4 


(20%) 


[19], [20] 


B[S -> Dri/]/iB[S -> L>£z/] 


0.416 ±0.128 






(31%) 


[22 E2] 


Am d 


0.507 ±0.005 


ps~ 


-1 


(1.0%) 


[16j 


Am s 


17.77 ±0.12 


ps~ 


-1 


(0.7%) 


[21 


r[Z 66]/r[Z hadrons] 


0.21629 ±0.00066 






(0.3%) 


[24J 


B[B -+X sl \/B[B -^X c eD\ 


3.346 ± 0.251 


io- 


-3 


(7.5%) 


L25J 



Table 1: Branching ratios used as inputs for the global 2HDM Type II analysis. They are listed and 
their values are given with their absolute uncertainty, their relative accuracy and the reference from 
where the value was taken. When two uncertainties are given, the first one is statistical and the second 
is systematic (often from theoretical origin). 



how strong the coupling of the meson to the axial current can be) and tm its lifetime. The 
corrective factor stands for channel-dependent electromagnetic radiative corrections. In 
this work, they are taken into account in the case of the lighter mesons (tt and K), where their 
impact is estimated to be at the level of 2 — 3% [2U US], and for the D meson, where the effect 
is 1% [SUES]- As far as -B-related observables are concerned, the experiments take into account 
soft photons corrections derived from their Monte Carlo simulated data, and we will assume 
that no further correction is required for these branching ratio (setting Sgf} = 0). 

The experimental accuracies for the branching ratios are given in the Table (TJ lying within 
~ 0.2—31% depending on the leptonic decay of interest. The main theoretical uncertainty arises 
from the decay constant, which is a non-perturbative quantity to be estimated by theoretical 
methods, such as quark models, sum rules, or lattice QCD (LQCD) simulations. We opt for the 
latter, since they provide well-established methods to compute these observables with a good 
accuracy and a satisfactory theoretical control. 

Over the last few years, many new estimates of the decay constants have been issued by 
different lattice collaborations, with different ways to address the errors. A part of the un- 
certainties has a clear statistical interpretation: lattice simulations evaluate correlators in an 
Euclidean metric expressed as path integrals using Monte Carlo methods, whose accuracy de- 
pends crucially on the size of the sample of gauge configurations used for the computation. But 
systematics are also present and depend on the strategies of computation chosen by competing 
lattice collaborations: discretisation methods used to describe gauge fields and fermions on 
a lattice, parameters of the simulations, such as the size of the (finite) volumes and lattice 
spacings used for simulations, the masses of the quarks that can be simulated, and the number 
of dynamical flavours (2 and 2±1 being the most frequent). In relation with these choices, 
the extrapolation of the results to physical parameters can be subject to different theoretical 
treatments (chiral perturbation theory, heavy-quark expansion. . . ), going beyond a naive linear 
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Input 


Value Unit Accuracy Reference 


Decay Constants 


Ik/ fn 
IdJ fo d 

fBsl ' fB d 

}b s 


1.205 ±0.0012 ±0.0095 (0.1%, 0.8%) [27J 
1.186 ±0.0048 ±0.0010 (0.4%, 0.1%) [27J 
246.3 ± 1.2 ±5.3 MeV (0.5%, 2.2%) [27] 
1.199 ±0.008 ±0.023 (0.7%, 1.9%) [27] 
228±3±17 MeV (1.3%, 7.5%) [27J 


Semileptonic Form Factors 


P 2 
A 

M v 

Mo) 


1.19 ±0.04 ±0.04 (3.3%, 3.3%) [25J 
0.46±0±0.01 (0,2.2%) [28] 
878 ±6 MeV (0.6%) [IS] 
0.9653 ± 0.0028 ± 0.0048 (0.3%, 0.5%) [27J 


BB mixing 


B Bs /B Bi 
Vb 


1.28 ±0.02 ±0.03 (1.6%, 2.3%) [27] 
1.05 ±0.01 ±0.03 (1.0%, 2.9%) [27] 
0.5510±0± 0.0022 (0,0.4%) [291 E0] 


Z ^bb 




0.02758 ± 0.00035 (1.3%) [24J 


b — > S7 parameterization 


C 

pole 

a s {m z ) 


0.546 ± ± 0.033 (0, 6.0%) [32J 
172.4 ±1.2 GeV (1.2%) [2S] 
0.1176 ±0.0020 (1.7%) [16J 


Running Quark Masses 


m u {2 GeV) 
m d (2 GeV) 
m s (2 GeV) 
m c (m c ) 
m b (m b ) 


2.40±0±0.90 MeV (0,38%) 16 
4.75 ±0± 1.25 MeV (0,26%) [16] 
96±0±30 MeV (0,31%) [IS] 
1.286 ±0.013 ±0.040 GeV (1.0%, 3.1%) [27J 
4.243 ± ± 0.043 GeV (0,1.0%) [25J 



Table 2: Parameters used as inputs for the global 2HDM Type II analysis. The parameters entering the 
calculations are listed and their values are given with their absolute uncertainty, their relative accuracy 
and the reference from where the value was taken. When two uncertainties are given, the first one is 
statistical and the second is systematics, often from theoretical error, hence treated in the Hfit scheme. 
For the latter systematics, whenever individual contributions are listed in the quoted reference we 
combine them linearly instead of quadratically, to stay consistent with the Rfit scheme. Therefore, 
our systematics errors can be larger than the one given in the corresponding references. Further note 
that the scale invariant top quark massrfit{mt) is computed from mP° le , following eq. (33) of 133^ with 
nf = 5 active flavours. 
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Figure 3: Leptonic decay of a meson through the exchange of a W boson (left) and a charged Higgs 
(right). 

extrapolation. 

The combination of lattice values with different approaches to address the error budget is a 
critical point of most global analyses of the flavour physics data, even though the concept of the 
theoretical uncertainty for such quantities is ill-defined (and hence is the combination of them). 
The CKMfitter group has collected the relevant LQCD calculations of the decay constants fs d , 
Ib 3 , Id s , Id, fx, fn (as well as bag factors Bb a , Bb s and Bk) and designed an averaging method 
aiming at providing a systematic, reproducible and to some extent conservative scheme [27J. 
These lattice averages are the input parameters used in the fits presented in this paper. 

In the specific case of light mesons (kaons and pions), the ratio of the decay constants fx/In 
is significantly better determined than the individual decay constants. It is hence worth con- 
sidering the ratio K^/^n of the kaon and pion leptonic partial widths instead of the individual 
branching ratios. It explicitly writes in the SM as [T5] : 

Y[K -> fiu] S M _ mK_ ( 1 - mf/m 2 K \ 
T[n ->■ /iz/]sm \1 -mf/ml J 

Let us notice that short-distance radiative corrections cancel out in the ratio. The long-distance 
corrections are accounted for by the parameter S^^ ni2 = —0.0070 ± 0.0035, which can be 
computed using chiral perturbation theory [31] . Similarly we also consider the ratio of taus to 
kaons to taus to pions decays. The latter writes: 

(f) 2 a+c 2/ n (6) 

The radiative correction is taken from [52] and reads 2 = 0.0003. The uncertainty is 

coming from 5^® 2 and ^^ 2 / Ke2 y ( - T7T2 / 7Te2 ) _ g_ggyg ^ 0.0027, both treated as iffit uncertain- 
ties. 

On the experimental side, the latest measurements of B — > tv branching ratio from the 
B factories [19, 20J have been taken into account. The two experiments BaBar and Belle find 
results in fairly good agreement and the weighted average of B[B — >■ tv] exhibits a departure 
from the SM prediction from the global CKM fit (more than 2a), due to tension between sin(2/3) 
and B[B — > tv], as pointed out in ref. [H]. 

In two Higgs doublet models, purely leptonic decays receive an additional contribution from 
charged Higgs, as illustrated in Figure [3j It turns out that in this case, this correction can be 
factorized from the SM prediction P [HI 36J: 

B[M -> tv] = B[M -> ev] SM (l + r H f, (7) 



1 1 (! + <W )• ( 5 ) 



Vud V /tt 



r[r -)■ Kv] SM 
T[r nv] SM 



1 — m l K lm l T 
1 — ml /ml 



ud 
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where the corrective factor r# writes in 2HDM Type II: 

r H = ( mqu - mqd tan2 P \ ( mM V (8) 
H V m Qu + m qd ) \m H +) 

A comment is in order concerning the structure of eq. (jHJ). Let us suppose that we have a 
perfect agreement between the measurement of a purely leptonic decay and its SM expression. 
There are actually in this case two distinct solutions in the 2HDM Type II framework: 

— th = 0, which can be obtained easily by sending m#+ to infinity. This decoupling solution 
corresponds to a general way of recovering SM predictions by assuming that all additional 
particles are very massive. 

— rjj = —2, which corresponds to a linear correlation between m#+ and tan/3. This fine- 
tuned solution depends on the mass of the meson and those of its valence quarks, and it 
is thus different from one meson to another. 

More generally, if we have a good agreement between the SM prediction and the measure- 
ment (as is the case for most leptonic decays), the 2HDM Type II fit will favour two regions: 
one region of high Higgs mass (related to the decoupling solution) and one diagonal band (in 
relation with the fined-tuned solution) in the parameter space (m H +, tan/3). 



2.2 Semileptonic B decay: B — > Drv 

Purely leptonic decays of mesons intertwine electroweak and strong interactions. However, the 
role of strong interaction boils down to the presence of a decay constant, to be assessed through 
theoretical methods. Semileptonic decays are more complicate to describe, since they involve 
form factors with a non-trivial dependence on the transfer momentum. 

If the form factors are known with a sufficient accuracy, semileptonic branching ratios start 
becoming valuable constraints on New Physics models - for instance, the comparison between 
leptonic and semileptonic decays provides a good test of the V—A structure of weak interactions. 
The BaBar and Belle experiments recently published first measurements of BiB — > Dtv) 

[2UI22]. 

An interesting observable is the normalized branching ratio TZb^Dtu — B[B — > Dtu]/B[B — > 
Deu] , which corresponds to a b — > c transition, with a CKM factor much larger than the purely 
leptonic B decay (and thus easier to study experimentally). In principle, the relevant form 
factors can be studied using LQCD simulations, as long as one is interested in a limited region 
of space-like transfer momentum, where both incoming and outgoing states consist of a single 
meson and no final state-interaction occurs. In order to extend the range of determination of 
form factors, fits combining lattice information and B — > Div can be performed [28], in order 
to constrain the shape of the relevant vector and scalar form factors over the whole kinematic 
regime. 

In the case of 2HDM-II models, the scalar form factor is a key ingredient, since it encodes 
the impact of the charged Higgs exchange in the semileptonic decay. Unfortunately, scalar 
form factors are notoriously difficult to handle on the lattice, and require dedicated methods 
and significant computing power to be estimated [31]. Due to helicity suppression, this scalar 
contribution arises with a factor mj/m 2 B in the amplitude, which means that only B — > Dtv t 
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is sensitive to this contribution 0. Following ref. p8] , we write the ratio IZb^dtv as a second 
order polynomial in the charged Higgs coupling, as: 



Kb^Dtv = a + a{JZe[s H ] + a 2 \s H 



(9) 



The factor sh accounting for the charged Higgs coupling was taken as: 



tan/3 2 




(10) 



sh = 



1 - m c /m b 



The polynomial coefficients ctj in eq. ([9]) depend on the slope parameter of the vector form 
factor p 2 and the scalar form factor A (taken as a constant) describing the B — ^ D vector and 
scalar form factors. A detailed statistical treatment of the theoretical uncertainties a la iffit 
requires an explicit knowledge of these dependencies. Therefore we do not use the final result of 
ref. [28J but we rather integrate their expression for the B — > Dtv differential branching ratios 
for various values of the parameters governing the form factors. We checked that the variations 
of the branching ratio with p 2 and A are smooth over the range of uncertainty of these latter 
parameters. Consequently, the coefficients cij are well parameterized as: 



where we have averaged over B& u — > D^u modes and where dp 2 = p 2 — 1.19 and dA = A — 0.46 
are the variations of the p 2 and A parameters. Let us mention that we have also investigated 
the constraints derived in ref. [38J from the same observables, and we checked that we obtained 
very comparable results, even though the theoretical expressions are quite different. 

2.3 Semileptonic kaon decay: K — > n£v 

The process K — > tt£i> is known to very good experimental and theoretical accuracies (better 
than 1%) and hence is included in this analysis, even though it involves lighter fermions. The 
branching ratio B[K° — > ^pv\j 'B[K° — > neu] writes: 



where 8 E '^ — (5.7 ± 1.5) • KT 3 and 5^ = (8.0 ± 1.5) • lO" 3 are radiative electromagnetic 
corrections estimated in Chiral Perturbation Theory in ref. [39J and recalled in ref. [15]. The 
phase space integrals I l K depend on both the scalar and vector form factors describing the 
K — > 7r transitions. In the low energy region of interest here, the vector form factor can 
be described accurately through resonance saturation (involving the K* pole). The scalar 
form factor is more delicate to describe, but it can be expressed through a dispersion relation 
involving data on ttK scattering. 

In principle, a similar analysis could be applied for B — > D*tv. However it involves four form factors which 
are poorly known and out of which only one would be sensitive to charged Higgs exchange |37j . 



a = 0.2970 + 0.1286 
ai = 0.1065 + 0.0546 
a 2 = 0.0178 + 0.0010 



dp 2 + 0.7379 
dp 2 + 0.4631 
dp 2 + 0.0077 



dA, 
dA, 
dA, 




(12) 
(13) 




(14) 
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W + H+ 

Figure 4: An example of b — >■ 57 transition through W exchange (left) and charged Higgs exchange 
(right). A second diagram can be drawn with the photon emitted from the charged boson. 

Exploiting these form factors, a parametrization of the ratio of the phase space integrals is 
derived in Appendix A: 

5 = Pl [z}(l+ P s[z}e G +p v [z}AM v ), (15) 
Ie K 

pi[z] = 7.96407 • 10 _ V + 2.59205 • 10" V + 7.82087 • I0~ 2 z + 0.647932, (16) 

p s [z] = 6.04408 • 10~V + 2.30011 • I0~ 4 z + 4.67096 • 10~ 4 z, (17) 

p v [z) = 4.24716 ■ I0~ 5 z - 2.69983 • HT 5 , (18) 

where AMy = (My — 892) is a parameter related to the mass of the K* resonance expressed 
in MeV. The factor ec is a free parameter in [—1; 1] reflecting theoretical uncertainties on the 
scalar form factor. Charged Higgs contributions are accounted for through the parameter z 
as PJ: 



f_K 1 | m s tan 2 g - m u / m 2 R0 - m\ + ^ 

f n f+(0) ° T m s -m u V m H+ 



where /+(0) denotes the normalization of the vector form factor at q = GeV 2 and Act £ 
[—0.5; 0] • 10~ 3 , estimated in Chiral Perturbation Theory, describes the deviation of the scalar 
form factor from fx/ fn at the Callan-Treiman point. 



2.4 Radiative decays: b — > sj 

The FCNC decay b — > 57 proceeding through penguin diagrams is a powerful benchmark 
to constrain the charged Higgs sector of New Physics models. The calculation of the B — > 
X s ~f branching ratio has been completed up to Next-to-Next Leading Order (NNLO) [JT] (see 
refs. 0Q1-021 for details). 

Its starting point is the effective Hamiltonian derived by integrating out the degrees of 
freedom heavier than the b quark, and expressed as products of effective operators (describing 
long-distance effects) with Wilson coefficients (including short-distance effects). The inclusive 
decay rate is expressed as the imaginary part of the relevant correlator of two 6 — )■ 37 currents, 
which can be expanded in powers of l/m&. Most of the effort in the field has been devoted to 
the computation of the perturbative part of the leading contribution to this expansion. 

The NNLO expression for the branching ratio is complicated; for practical purposes, we have 
chosen to parametrize the results given by the public package SusyBSG [13] (other programs 
do exist, see for instance [S]), based on Leading Logarithm (LL) expressions. The SusyBSG 
code includes NLO perturbative corrections for different theoretical frameworks: SM, 2HDM 
and MSSM. 
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In the 2HDM Type II, the exchange of charged Higgs bosons in addition to charged weak 
gauge bosons, shown in Figure HJ provides a further contribution to the relevant Wilson coeffi- 
cients of the effective Hamiltonian. The parameters of the SM prediction, limited to NLO, have 
been tuned in order to recover the most accurate NNLO result given in ref. [41] . Following the 
notation in ref. [45J, the normalized branching ratio for B — > X s j writes: 



_ B[B X sl ] 



v t *v tb 



V, 



do 



^(P + N), (20) 



B[B ->■ X c tv\ 

where P denotes the leading contribution in the l/m& expansion, computed in perturbative 
QCD, and N the non-perturbative ones, corresponding to higher orders in the l/m& expansion 
(starting at 1/m 2 ). If P can be systematically improved by computing higher and higher orders 
in perturbation theory, it is hard to provide more than an order of magnitude for N [4"U] . 

The normalization factor C accounts for the phase space difference between charmed semilep- 
tonic transition and B — > Xgj decay. For our analysis of 2HDM Type II, P and N are 
parametrized using two functions, A and B, depending on a reduced set of the relevant input 
parameters (a s (mz), Tnf ole and m c (m c )). Making use of the perturbative expression of P at 
the leading-logarithm level, the functions A and B are defined as: 

P + N = (C$g> + BAC; S H { ?) 2 + A, (21) 

and fitted to reproduce the results from the SusyBSG package. In eq. (l2~Tj) . the factor AC^+ 
models the Charged Higgs contributions. A and B, which are independent of the 2HDM Type 
II parameters and tan/3, exhibit smooth linear variations with the input parameters. 

Further details on the parametrization used in this analysis as well as the formulae for A and 
B functions are given in Appendix B. 



2.5 Neutral 5-meson mixing 

In the Standard Model (SM), neutral meson mixing occurs due to box diagrams with two W 
exchanges. In the case of Bd and B s mesons, the hierarchical structure of the CKM matrix and 
the large mass of the top means that the mixing is dominated by short-distance physics coming 
from diagrams where the internal fermion lines are top quarks. In two-Higgs doublet models, 
the observables related to neutral-meson mixing receive charged Higgs contributions [31 H] 151 152"] . 
Indeed, one gets further diagrams obtained by replacing one or two W lines by a charged Higgs, 
yielding 



G 2 

Am, = -^^{V tq Vt b ) 2 ri B m B m 2 fl q B Bq (S ww + S WH + S H h), (22) 

Q fu 9 6 6x 2 tw \n(x tw ) \ 

V 1-XtW [l-Xtw) {i-x tw y J 



x t H f (2x w - 8) ln(x tH ) 6x HW \n(x tw ) 8 - 2x tw 

tan 2 f3 \(1 - x HW )(l - xm) 2 (1 - x H w)(l - Xtw) 2 (I - x tw )(l - x t H, 



x tH ( l + x tH 2x tH \n(x tH )\ 



tan 4 /3 \(l-x tH ) 2 (l-x tH y J 

with f/m 2 . Sy/Wi Swh and Shh indicate respectively the internal bosonic lines of the 

corresponding diagrams with an external light quark q = d,s. Analyses including radiative 
corrections are available [17], but the leading-order expressions above are sufficient for the 
required accuracy in our present purposes. 
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2.6 Constraint from electroweak precision data: the Z — > bb vertex 

It is time to make an excursion at the border of flavour physics. The Z bb vertex has provided 
opportunities to search for New Physics contributions, due to the heavy masses involved. In 
particular, the radiative corrections at the vertex might imply charged Higgs exchanges in 
addition to the standard Wtb couplings. The partial width T[Z — > bb} is subject to sizeable 
QCD corrections. It is hence relevant to define the ratio of the Z partial widths Rb = T[Z — > 
bb]/T[Z — > hadrons] for which most QCD corrections are suppressed. Rb has been measured 
by the LEP experiments with a remarkable accuracy [21]. We will not consider neutral Higgs 
corrections here. Indeed it has been shown [IB] that they contribute significantly for large values 
of tan /3 only, where Rb is not a competitive observable with respect to the other observables 
considered here. On the other hand, the Rb measurement can yield valuable constraints to 
exclude regions of the 2HDM Type II parameter space at low tan /3. 

Following the work described in refs. [HI H9J [50] , we parameterize Rb as: 

1 K b 

R b " 1 + (gt - W + + W 

L Ft 

where g h ' are the couplings of left- and right-chirality b quark to the Z boson; Kb is the sum 
of the axial and axial-vector couplings of the quark flavours lighter than the b quark, embed- 
ing QCD and QED radiative corrections remaining in Rb as well as effects from the 6-quark 
mass. The 3 latter SM quantities have been parameterized to reproduce the predictions of the 
ZFITTER package [51] . which depends primarily on the top quark mass, m t . However, one 
has to take care about the correlations between the top quark mass and the rest of the elec- 
troweak parameters. In particular, the neutral Higgs mass, m H o, is constrained by a dedicated 
Electroweak fit [23], following High-Q 2 data fit, but excluding the direct measurements of Rb, 
A°p b B and Ab from the fit inputs. The dependency of the Electroweak % 2 on mt is modelled 
by a correlation factor between m t and m#o. The full details of the parametrization are pro- 
vided in Appendix C. With the top quark mass quoted in table [21 it yields the SM prediction 
Rb = 0.21580(4), excluding the direct measurement of Rb as well as Ap b B and A b from the fitd. 

The charged Higgs contribution induces a redefinition of the coupling constants (ref. [IB] 
corrected in [50]) according to: 

-l _ -L | G F m 2 w f m t 1 V 

9b ~~ 9b,SM t y- p J " z 

8V27T 2 \m w tan/} J 
- n R_- n R °F m w ( m b ,n\ 2 F 

9b - 9b,SM ~ " nz 9 -— tan " tz 

8V27T 2 \m w ) 

The function F z takes into account two-loops corrections following the work developped in [SD] 
and a parameterization is provided in Appendix C. 

Following [19], we assume that the oblique corrections due to the second Higgs doublet are 
negligible, so that the modifications from the 2HDM will mainly affect the vertex corrections, 
and thus modify b-quark observables (and hence T(Z — > bb)) before any other observable. We 
are then still allowed to determine the SM prediction for Rb from the Electroweak fit described 
above (which does not include any b- related observables) , out of which we can deduce the value 
of Rb in 2HDM using eqs. [23 

2 Let us note incidentally that this prediction is in excellent agreement with the ZPOLE Fit results of [24] . 
which include the measurements of Rb, A°p b B and A^. This is due to the fact that the direct measurement of mt 
coincides with its prefered value from ZPOLE only inputs, which significantly depends on Rb. 



m: 



m 



(27) 



mt 



in 



H 4 



(28) 
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3 Individual constraints in the (ra#±,tan/3) parameter 
space 

All the fits reported in this paper are performed within the frequentist statistical framework 
advertised in ref. [10]. In all analyses, the CKM matrix parameters are determined simultane- 
ously with the 2HDM additional parameters. Hence, we recall first how the CKM parameters 
are measured in the framework of the SM and how their determination is modified once the 
2HDM hypothesis is tested. 



3.1 Standard Model inputs and parameters 

There are four free parameters of interest describing the CKM matrix (A, A, p and ff) in the 
extended Wolfenstein parametrization [TO] : 

A= ; Ksl =, A\*= lVcbl =, p + i fj = -Wjk. (29) 

V ud \ 2 + \V US \ 2 \J\V ud \ 2 + \V us \ 2 cdcb 

The dependence of the other CKM matrix elements on these parameters follows from the 
unitarity of the CKM matrix. This definition ensures unitarity at all order of the development 
in power of A and warrants the p and fj parameters not to depend on phase conventions. 

In the SM, A and A are accurately determined: A is measured from super-allowed nu- 
clear transitions and semileptonic kaon decays and A comes from the inclusive and exclusive 
semileptonic 6-hadron decays with charm. On the other hand, the parameters p and fj, being 
respectively the real and imaginary coordinates of the unitarity triangle (UT) apex, are less 
constrained. The fit of the CKM matrix in the SM hypothesis and the metrology of its four 
parameters make use of several observables: 

- \V u d\, \V US \ and |Vrf,| determine the A and A parameters and fix accordingly the length 
scale of the UT. 

- \V u b\ (including B(B — > tu)), Am^ and Am s are CP-conserving observables, sensitive to 
the sides of the UT. 

- a, 7, sin 2/3, cos 2/3 are CP-violating observables measuring the UT angles from 5-meson 
decays whereas |e#| assesses CP violation in kaon mixing. 



3.2 2HDM inputs and parameters 

Let us move to testing the 2HDM Type II hypothesis. It requires to add the two 2HDM Type 
II parameters (m H +, tan/3), but also to modify the set of constraints used in the global SM fit 
to fix the CKM matrix. We have therefore to split our observables into those used to fix the 
CKM matrix, and those needed to constrain the additional parameters from the 2HDM Type 
II. 

First of all, the observables dealing with neutral-meson mixing proceeding through AF = 2 
transitions will receive charged Higgs contributions [3J HJ El E2], and we can use neither the 
oscillation frequencies of B d and B s mesons (Am^ and Am s respectively) nor the CP-violating 
parameter e#, as inputs for the CKM matrix. Moreover, the UT angles a and cannot 
be used independently, since their determination relies on an interference between decay and 
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mixing. However, it turns out that the combination of a and /3 inputs in the fit constitutes a 
determination of the angle 7 in which AF = 2 contributions cancel [53J. 

Moreover, the AF = 1 processes proceed through a jy-exchange in the SM, and they are 
thus affected directly by the presence of a charged Higgs boson. However, we know that these 
contributions are proportional to the masses of the quarks and charged leptons involved. At 
low energies, the exchange of a charged Higgs boson yields four-quark operators with weights 
m i ' m 2/Mfl, where mi and 7712 are the masses of two fermions (quarks or charged leptons) 
involved in the four-quark operator but not coupled together with a Higgs, as can be seen from 
eq. ©. 

We expect therefore that only processes involving massive quarks and leptons will be very 
sensitive to 2HDM Type II contributions, which selects naturally some of the processes con- 
sidered above for the determination of (m#+, tan/3): p or t leptonic decays of B, D s and D, 
B — > D p or r semileptonic decays (tree processes), Z — > bb, b — > 57 and neutral B meson 
mixing (processes with a top-quark loop). 

On the other hand, the CKM matrix, and thus the apex of the unitarity triangle, can be 
determined by taking AF = 1 processes where at most one heavy mass is present. This selects: 

— the determination of 7 from a + /3 (rrib ■ m u /m 2 H , ■ md/m"^), 

— the determination of \V cb \ from semileptonic b — >■ c decays (m b ■ m e /m 2 H , m c ■ m e 

— the determination of \V u b\ from semileptonic b — > u decays (m?, • m e /m 2 H ), 

— the determination of \ V u d\ from super- allowed /3 decays of nuclei (no heavy mass involved). 

The determination of 7 from B — > DK does not enter this list as it scales like nif, ■ rris/Mfj. 

Figure [5] shows the combined constraint in the (p, ff) parameter space. Though less con- 
straining than the SM global fit, the constraints chosen here yield two well-defined symmetrical 
solutions for the apex of the unitarity triangle. The achieved accuracy is mostly due to the 
world average a determination, driven by the latest B — >■ pp measurements (see ref. [27] and 
references therein). 

Before discussing the combined analysis of all the above constraints, we would like first to 
focus on the most stringent individual constraints in the parameter space (m#+,tan/3) among 
the different classes of observables potentially sensitive to charged-Higgs exchanges. If we 
impose the Higgs sector to remain in a perturbative regime, an upper bound on the value of 
tan/3 can be obtained around 200 [5], and our plots will correspond to this region (with a 
logarithmic scale). Two different constraints turn out to constrain the 2HDM Type II very 
efficiently: the leptonic decays and the b — > 57 branching ratio. 

3.3 Leptonic and semileptonic decays 

As mentioned in the introduction, the only observable which departs significantly from the SM 
prediction is B(B + — > t + v). According to eq. ([8]), 2HDM Type II can accommodate such a 
sizeably larger value with respect to the SM in only two situations: 

3 We should rigorously have taken solely the electronic semileptonic b —> c (and conversely b — > it) decays. 
Yet, the electron and muon averages we are producing are by far dominated by theoretical uncertainties. On 
top of that, electronic and muonic extractions are in very good agreement and the split averages are not easily 
available. 
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Figure 5: Superimposed individual constraints for the fit comprising the observables which involve 
light fermions at 95% CL (and excluding AF = 2 observables). The yellow area is the solution driven 
by the combination of individual constraints at 95% CL. The unitarity triangle drawn here is obtained 
from the global fit displayed in Figure^ 



— at low enough tan j3 with respect to y^m u /rrib — 0.02, to get a positive correction ru to 
the SM value, 

— in a fine-tuned scenario (defined in section 2.1) with r# < —2, so that the 2HDM Type 
II contribution (1 + r#) 2 enhances significantly the SM prediction. 

The 95% CL constraints derived in the plane (m#+, tan/3) from the various leptonic decays 
are shown in Figure [6j in a log- log scale. Let us recall that the charged Higgs contributions are 
identical for K^/it^a and r — > Kv/t — > ttv. However, since the experimental value of K^/ff^ 
is yet well better known than r — > Kv/t — > ttv, the latter only has a marginal contribution to 
the combination. 

Combining all constraints from leptonic decays, the minimum \ 2 value of Xmin = 14.8 (p- 
value = 85.5 ± 0.3%) is found at small where B — > tv and D — > lv fine-tuned regions 
overlap. The very small charged Higgs mass, excluded by direct searches at LEP, reflects the 
fact that 2HDM Type II can hardly accommodate the large value of the measured B — > tv 
branching ratio at low masses of the charged Higgs but in a fine-tuned scenario, as can be seen 
in Figure [71 

At values of tan f3 sufficiently large compared to a/ m qu / m qi , the H ± contribution to the 
branching ratio behaves as r# ~ —(tan 2 (3)m* s /m 2 H+ with m c g = fflM\/%/(% + m g d )- Since 
the corrections to the SM branching ratio are at least quadratic in the effective mass m e fj, the 
B — >■ tv branching ratio sets a constraint on the high-m#+ region. Let us mention that the ratio 
K^/t^, accurately measured and predicted, plays a significant role in the combined leptonic 
limits. As shown in Figure UJ this additional observable disfavours the B — > tv fine-tuned band 
in the region tan (3 > 10. 
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Figure 6: Constraints on 2HDM parameter space (m H +, tan/3,) from purely leptonic decays. The 
upper left plot stands for constraints from K — >• \iv and r — > Kv decays, the upper right for D — >• fiv 
decays, the lower left for D s decays and the lower right for — > rv. The superimposed green line 
delimits the la confidence area. The excluded regions (white parts of the plots) correspond to more 
than 95%CL. 
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Log iQ [ tan(P) ] Log J tan(P) ] 

Figure 7: Combined constraints on 2HDM Type II parameter space {m H + , tan (3) from purely leptonic 
decays of K, D and B mesons. The different sets of observables are defined in the left Figure. Darker 
to lighter shades of colours correspond to (Knil^ti +t — > Kv), (D — > \iv), (B — >• tv) and (D s — > iv) 
individual constraints. The complementary area of the colored region is excluded at 95% CL. For the 
sake of clarity, the figure on the right displays the combined constraint alone at 95% CL. The dark 
green line defines the la confidence area. 
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Figure 8: Combined constraints on the ratio m H + / tan f3 derived from leptonic and semileptonic decays 
of K , D and B mesons, in the large tan j3 limit. 
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The limits derived from the semileptonic decays B — > Drv and K — > n£v are similar in 
shape to those derived from the purely leptonic decay B — > tv, although less constraining than 
the latter. Hence there are not further displayed. Note that despite BiB — > Dtv) suffers from 
larger theoretical and experimental uncertainties, because of the heavier mesons involved in the 
decay, it is 2 times more constraining on 2HDM(II) parameters than K — > 7r£v. 

It is very illustrative to compare all constraints derived from leptonic and semileptonic 
decays in the large tan (3 limit (tan (3 > 30). In this limit, the charged Higgs contributions depend 
on a single coupling parameter, namely tan f3 /m H +. Figure [S] shows the limits on m H + / tan/3 
derived from the most challenging observables, i.e., B — > tv, K^/n^ and B — > Dtv. We also 
show the combined limit from all leptonic and semileptonic decay observables considered in 
this study. Both and B — > Drv turn out to exclude the fine-tuned solution due to 

B — > tv, at more than 95% CL. It results in the combined limit m#+/ tan/3 > 13.1 GeV at 
95% CL, from all semileptonic and leptonic decays in the large tan/3 limit. 

3.4 Loop processes 




Log iQ [ tan(|3) ] Log J tan((3) ] 

Figure 9: Superimposed constraints on 2H DM parameter space (m#+, tan/3j from BB mixing and 
Z — » bb (on the left) and from from b — > S7 branching ratio. The colored area are confidence regions 
at 95% CL. 

The constraints derived from BB mixing and Z — >■ bb are similar in shape as can be seen on 
Figure [9j The constraints derived from Am^ and Am s are twice as stringent as those arising 
from Z — > bb. They all exhibit a divergence with 1 / (m^+ tan /3) resulting in small values of 
tan (3 to be disfavoured except at very large values of m#+ . 

The b — > X s 7 branching ratio, where X s denotes any charmless final state with strangeness, 
is measured using either semi-inclusive or inclusive method. The Heavy Flavour Averaging 
Group (HFAG) proposes an average of measurements (2S], which is used in the present work, 
performed by the CLEO [54J, Belle [55] and BaBar [56] experiments (only the reference of 
the latest measurements is given for each collaboration). The key point of the branching 
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ratio determination is the correction implied by the photon energy threshold experimental cut 
and derived from a theoretical model of the photon spectrum shape. In that respect, HFAG 
advocates the use of the extrapolation factors determined in [57]. The average is: 

B{b -> sj) = (3.52 ± 0.23(stat) ± 0.09 (syst)) 10" 4 . (30) 

The b — > 57 branching ratio is currently the constraint which dominates the global 2HDM 
Type II fit. Figure |9] shows the exclusion region at 95% CL in the parameter space (m#+, tan /3). 
This result is in fair agreement with other determinations [32]. Let us mention that two 
critical parameters in the determination of this limit are the charm quark mass fn c (m c ) and 
the semileptonic phase space factor C. As an illustration, if the central value and uncertainties 
of m c are varied from the one quoted in Table [2] to the one used in ref. [31] (fn c (m c ) = (1.224 ± 
0.017 ± 0.054) GeV), the limit on m#+ is increased by 10%. Similarly for the parameter C, 
if to use the central value and errors from ref. [31] (C = 0.58 ± ± 0.016) the limit on m#+ 
is decreased by 9%, despite the error on C is two times smaller than the one we used in our 
analysis. 



4 Combined Analysis 
4.1 Goodness of fit 

In order to test the relevance of the 2HDM Type II all observables are compared to their 
theoretical predictions within a combined global x 2 test. The \ 2 is minimized over all the- 
ory parameters, and theoretical uncertainties are treated in the iffit scheme as for the global 
CKM fit [10]. There are 76 observables (the determination of the a angle corresponds to 41 
observables alone) for 56 parameters, yielding xLnOtsf^HDM) = 22.54. From a Monte-Carlo 
toy-experiment study, the corresponding p-value is found to be p(2HDM) = (64.8 ± 0.8)%, 
assuming for the true parameter values those found during the initial x 2 minimisation (plug-in 
p- value scheme [ID]). It is worth comparing the observed Xmin value in the 2HDM scheme to 
the one obtained for the same observables but with their SM predictions. For the latter, the 
Xmin,obs * s Xmin,obs(SM) = 22.56 and a corresponding p-value p(SM) = (69.1 ± 0.8)% is found. 
Let us stress the Xmin vam e can only decrease when we move from the SM to 2HDM predic- 
tions: the observables are the same but the 2HDM Type II has additional parameters which 
can reproduce the SM predictions in the particular limit m#+ — > oo. Therefore, by comparing 
the observations to the predictions alone, one cannot reject or nullify the 2HDM Type II while 
keeping the SM. Nevertheless, the almost equal value of Xmin 101 both models leads us to the 
qualitative conclusion that the 2HDM Type II does not perform significantly better than the 
SM. 

We can measure more quantitatively how the agreement with the data improves once one 
moves from the SM to the 2HDM Type II. For this purpose, we introduce a new test statistics, 
AxLn, defined as: 



A Xmin = Xln(SM) - Xmm( 2HDM )" (31) 

High values of ^Xmin would indicate a deviation of the observation from the SM prediction 
that could be accommodated by 2HDM Type II. The cumulative distribution of this test is 
derived from a Monte-Carlo toy-experiment study where the SM parameter values have been 
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Figure 10: IS.x 2 min cumulative distribution within the SM. The black squares stand for the distribution 
obtained from Monte-Carlo assuming the SM as the truth. For comparison, the red dots show the 
distribution obtained for a 1 degree of freedom x 2 distribution with identical binning. 



fixed to the values found in the global minimisation (this is the truth hypothesis to nullify). 
The result of this toy analysis is shown in Figure [TUJ A departure of ^-Xmin at three standard 
deviations (a p-value of 0.27%) would correspond to A%^ in > 9.2. The observed value A%^ in = 
0.02 amounts to a p-value of ~ 100%. Therefore, the toy analysis does not give any reason to 
reject the SM in favour of the 2HDM Type II. Conversely, the 2HDM Type II hypothesis is not 
invalidated and the next section is dedicated to the derivation of the exclusion limits on the 
model parameters (mfl-+,tan/3). 

Eventually, we point out that the test distribution was found to be well approximated 
by a x 2 distribution with one degree of freedom while one would have naively expected two 
degrees of freedom brought by the two additional free parameters m H + and tan/3. This could 
be understood from the fact that the fit is dominated by the b — > 57 constraint, almost one- 
dimensional and depending only on m#+. 

4.2 Combined limits on 2HDM Type II parameters 

Figure [TT] shows the combined 95% CL confidence area in the plane tan (3). The minimum 

X 2 of Xmin — 22.5 is obtained for m H + ~ 4 TeV. At high H + masses and irrespective of the 
value of tan/3 (decoupling limit), the charged Higgs contribution becomes negligible for all the 
processes we are considering in this analysis, so that the SM predictions are recovered. 

At large tan/3 ( > 30), leptonic (mainly B — > tv) and semileptonic constraints compete with 
b — > 57 and sharpen its exclusion limit. At small tan/3 (^$1), the most stringent constraint 
arises from the BB mixing and to a second extent from Z — > bb. These results can be compared 
with those obtained from the Gfitter group [12J, which performed a global fit to electroweak 
precision data both in the Standard Model and in the 2HDM Type II. In the latter case, the 
observables involved were B — > iu, B — > Drv and a ratio involving kaon decays (namely K^, 
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Figure 11: Global constraints on 2HDM parameters m H + and tan/3 from all analysed observables. 
Each color corresponds to a different set of observables, as quoted in the Figure. The complementary 
area of the colored one is excluded at 95% CL. The horizontal black line indicates the 95% CL limit from 
direct searches at LEP }58\/ . The dotted line within the orange combined area delimits the corresponding 
la confidence area. For the combination of leptonic and semileptonic constraints (light green area) we 
assumed that the p-value is best approximated by a 1 d.o.f. \ 2 distribution since 2HDM contributions 
essentially depend on the ratio mjj+/ tan(/3). 
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Figure 12: One-dimension constraint on mjj+. 



K^2 and 717^), b — > 57, Z — )■ 66 (with a more dedicated treatment of the latter observable than 
in our work). The values of CKM matrix elements involved in flavour observables were taken as 
external inputs, whereas we determine them from the fit described in sec. 13.21 We notice rather 
similar exclusion areas for the various individual constraints (e.g., the existence of fine-tuned 
solutions for leptonic decays), apart from a slightly different shape in the case of Z — > bb. As 
in our case, b — > sj favours high values of the charged Higgs mass, irrespective of the value of 
tan f3. 

Figure d2] shows the one-dimension constraint found in the global analysis. A lower limit of 
the charged Higgs mass can be inferred: 

m H+ > 316 GeV at 95% CL [this work] 

while no significant constraint is obtained for tan /3. 

It is interesting to compare these results with the bound obtained from direct searches at 
LEP for any value of tan/3 [58] (specific studies are also reported in refs. [59]): 

m H+ > 78.6 GeV at 95% CL [direct] 

5 Conclusion 

In the last decade, the B factories have performed a set of remarkable measurements in the 
quark sector which yielded an impressive overall agreement between the SM and the data. 
Tight constraints on the theories beyond the SM can be inferred: flavour-changing neutral 
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currents are small as predicted in the SM and the KM mechanism has proven to describe 
the observed CP-violating phenomena in flavour physics with a good accuracy. Among the 
theoretical extensions of the SM, the addition of second scalar doublet (2HDM Type II) is 
particularly appealing, so that the flavour structure of the SM is preserved. 

In this article, we have discussed tests of the 2HDM Type II in the light of recent flavour 
physics data. In particular, the measurements of purely leptonic decays of B and D mesons 
(which departed more or less significantly from their SM predictions) have been combined to 
obtain a comprehensive combined constraint. They have been analyzed together with comple- 
mentary constraints such as semileptonic decays, b — >■ 57 or Rb measurements. 

The outcome of this combined analysis is that the 2HDM Type II is not favoured by low- 
energy data due to the interplay between B[B — > t + u] and b — > S7 measurements, or at least, 
that it does not perform better than the SM. If we assume that the 2HDM Type II is realized in 
Nature, constraints on its parameters can be derived from this global analysis. In that respect, 
a limit on the charged Higgs mass m#+ > 316 GeV at 95% CL is obtained irrespective of tan 0. 

This analysis considered tree and loop-induced AF = 1 processes, as well as AF = 2 mixing 
processes for neutral B mesons. Loop-induced rare decays such as B s — > fi + fi~ might receive 
contributions from an extended Higgs sector, and would provide further constraints on this 
model. A natural extension of this work would include such observables in order to perform a 
more comprehensive test of the Two Higgs Doublet Model of Type II. 



Acknowledgments 

We thank our colleagues of the CKMfitter group for many stimulating discussions and for a 
critical reading of the manuscript. We also thank J. Orloff and N. Mahmoudi for very useful 
comments on two Higgs doublet models. We eventually thank H. Flacher and A. Hocker (from 
the Gfitter group) for valuable discussions on the constraints and P. Slavich and G. Degrassi 
for their useful comments on the Rb observable description. This work was supported in part 
by the ANR contract ANR-06-JCJC-0056 and by EU Contract No. MRTN-CT-2006-035482, 
"FLAVIAnet" . 



A Parametrization for K° — > nlis 

The phase space integrals I l K in eq. ( 1T4|) is split into two integrals as l l K = i l K + + J^ : 



2 



Iko = f (1 - rl)n f"^ *M(1 - r }f /oM . (33) 



involving the rescaled transfer momentum u = t/m 2 K0 , the rescaled masses rv = (m^/m^o) 2 , 
r£ = [mf/m K o) 2 and the phase space term A [it] = a/(u — (1 + rv) 2 ) — 4r 2 . 

A representation of the two Kix form factors is needed to compute these integrals. The vector 
form factor / + is not sensitive to charged Higgs contributions. We use a pole parametrization 
with pole mass My tuned on experimental data. For the scalar form factor, which is sensitive 
to H + contributions we follow the dispersive parametrization given in ref . [60] . This is based on 
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a dispersion relation proposed with two subtraction points where the scalar form factor is well 
known: the normalization at (determined from lattice QCD) and the Callan-Treiman point 
(at t = m 2 K — m 2 ) where it is predicted to be related to fx/ f-n up to very small corrections. 
It requires the knowledge of ttK scattering phase shifts at low energies, obtained through the 
analysis of Roy-Steiner equations (ref. [6l]). 



i+M = ry , M u ] = ex P ( i - [z-g 

r v — u \l-r 7 



u 



(34) 



with the rescaled mass of the vector resonance ry = (Mv/rriKo) 2 , and the subtraction constant 
of the dispersion relation 

f f If If fn\j_A ms tan2 & ~ m » ( m K° ~ m l+ ^ fo K \ 

z = fo = fx U //+ + Act 2 35 

m s -m u \ m 2 H+ J 

where charged Higgs contributions are included. 
The functional dependence of G is: 

G[x] ■ 10 2 = (2.09 + 0.26e G )x + (3.98 + 0.56e G )(l - x) + (0.45 + 0.01e G )ar(l - x), (36) 

where we introduced theoretical uncertainties through the parameter e G lying in [— 1;1]. 
Substituting the form factor expressions in eqs. ( 132|) and (|33|) . the ratios of phase space integrals 
are computed by numerical integration for various values of the parameters z G [—0.4; 0.8], 
e G G [—1; 1]. The relative variations of the integrals ratio propagate linearly with e G and AMy 
over the investigated range. The local derivative with respect to the latter parameters only 
depend on z and their variations are smooth. These observations led us to the parameterization 
given in eq. (ITS]) , in perfect agreement with our numerical studies. 



B Parametrization for b — > 57 
B.l Analytic expressions 

At the order of leading logarithms, the relevant Wilson coefficient for the effective SM Hamil- 
tonian, C^gffl, of eq. ( 12"T|) writes as: 

C e 7 !sV = V%Cl%\M + - V^lM + E M", (37) 



i=i 



where rj = a s \po]/a s \pb) is the ratio of the strong coupling constants at the scales /io and 
fib- We followed ref. (13] by setting fi = 2mw and /x& = 1.5 GeV as the mass scales for the 
renormalization scheme. The Wilson coefficients C^ SM and Cg ^ M , at scale fi t , write: 

= -^ L (2F 1 [x m ] + 3F 2 [x tw }), (38) 
<LW = -^fiM, (39) 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


hi 
di 


626126 

272277 
14 

2a 


56281 

51730 
16 

23 


3 

6 7 
23 


1 

14 
12 

23 


-0.6494 
0.4086 


-0.0380 
-0.42230 


-0.0185 
-0.8994 


-0.0057 
0.1456 



Table 3: Magic numbers required to evaluate the LL perturbative contributions to b —> s, taken 
from ref. fjhj . 



with x t w = (n^t[^o]/mw) 2 and the functions Ff 



x 3 — 6x 2 + 3x + 2 + 6x ln(x) 

FlW = m^TT* • < 40 > 



^ r . 2x 3 + 3x 2 — 6x + 1 — 6x 2 \n(x) 

Fzlx] = - — . (41) 

1 1 12(x- l) 4 V ; 

The running of a s is computed at the second order with rif = 5 active flavours. One has: 

a\ r a i 127r A 3481n(L A )\ 

with Z m a = 2 Kqcd) and the matching condition «s [mz, Aqcd] = ot s {mz). Moreover, the 
running quark mass, m t [//o], at scale /i is computed at 3 rd order in a s from the pole quark 
mass, following [33J. For nj = 5 active flavours, it goes as: 



mM = m t pole (l - ( - + ) a, - ( 9.125 + — Z„ M + -Zj M a\ (43) 



4 , \ / 419, 2 

3 + W ) <V - (9- 125 + 72" W + g' 

- (0.3125Zj M + 4.5937Z 2 M + 25.3188Z mM + 81.825) a 3 ) (44) 



with l^ M = 2 ln(/i/mf ole ) and a M = ai [/i, Aqcd] /ft- Finally, the third term in eq. ( 157)) further 
accounts for the remaining Wilson coefficients Cj = [i ; 6]. We collected these relevant numbers in 
Table El 

Charged Higgs contributions play in at leading logarithm order in the Wilson coefficients 
and Cg '' by substituting the W ± boson with the charged Higgs H ± . They result in the 
substitution in eq. ([37]) : C$ M -> Cf°j A/ + ACf^ + , with: 

AC^M = - X -f- ^-^^-F x [x m \ + F 2 [x tH \) + + F 4 [x tH ]^j , (45) 

where = (mi[/xo]/m^+) 2 and F 3j4 is defined as: 

_ x 2 -4x + 3 + 21n(x) _ x 2 -l-2xln(x) 

3[J " 2(a;-l) 3 ' 4[J " 2(x-l) 3 ' [ ' 

Following, the charged higgs contributions, AC^+ , to the effective Wilson coefficient fall 
from eq. (|37j) . as: 

AC^f = V %C*% + \M + ffo* - (48) 
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X 


x c 


x t 




X s 


A 


3.155 • 1(T 2 


2.80 


-1.06 • 


10~ 4 


36.2 


B 


7.564 ■ KT 1 


-2.43 ■ lO^ 1 


-7.68 ■ 


10" 4 


-4.62 



Table 4: Coefficients of the Taylor expansion of A and B functions. The label X stands for A or B. 



The higher-order perturbative terms contain charm quark contributions and electroweak 
effects in addition. Hence, they also depend on the running quark mass m c (m c ) and on the 
photon energy experimental cut-off, Eq ~ 1.6 GeV. They are known exactly at NLO and 
parameterized at NNLO [41] . A full expression for the NLO term is very lengthy, as can be 
inferred starting from the final results given in ref . |JT] • The non perturbative corrections seem 
to be only partly known [H)]-[II]. They are small as compared to other terms (a 10% effect) 
and uncertainties mainly depend on m c |45j . 

B. 2 Parametric model 

Perturbative (P) and non-perturbative (N) contributions are modeled according to eq. (T2"Tj) 
with two parametric functions A and B, depending on a subset of input parameters: m^ ole , 
m c {m c ) and a s {mz)- We neglect the dependence on the photon cut-off energy, which is not 
quoted in ref. |41] . Following our assumptions, the values of A and B were computed for 
different sets of parameters by fitting the B — > X s 7 branching ratio given by SusyBSG as 
a function of the charged Higgs mass, m#+. We found perfect agreement between SusyBSG 
results and our model, which justifies a posteriori the expression given in eq. (1211) and our 
parametric assumptions. 

We further studied the variations of A and B as a function of the inputs m^ ole , m c (m c ) 
and a s (mz)- We cross-checked that the variations of A and B with respect to the others 
parameters, quoted in Table 4 of ref. [H] , are negligible. Those with respect to m^ ole , m c (m c ) 
and a s (mz) are linear for small changes. We notice that linearity holds for deviations of the 
input parameters as strong as 3a from their central value. Therefore, it seems very appropriate 
to approximate the functions A and B by their leading-order Taylor expansions around the 
central input parameter values: 

A = A (l + a c Am c (m c ) + atAm^ ole + a s Aa s (mz)), (49) 
B = B (l + b c Am c (m c ) + b t Am p t olc + b s Aa s (m z )), (50) 

where Ap denotes the variation of the input parameter p around its central value. The coeffi- 
cients values were extracted from a polynomial fit, summarized in Table SJ 

C Parametrization for Z — »• bb 

C. l SM prediction 

Following [21], we performed a global Electroweak fit to the Z pole observables in the High-Q 2 
scheme, assuming purely Gaussian errors and taking their correlations into account. There are 
five fit parameters: the masses of the Z boson, the neutral Higgs boson and the top quark as well 
as the strong coupling constant, a s [mz], and the hadronic vacuum polarization, Aa^ d [mz], 
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considered at the Z pole. The values of the Z pole observables and the hadronic vacuum 
polarization are taken from [21] (table 8.4) but excluding the measurements of Rb, A°p b B and At,. 

The top quark mass is essentially constrained by the direct measurement given in Table [2] 
while the W boson mass and width are taken from JTB] . With these inputs we get: Rb = 
0.21580(4) and log 10 (m H o/GeV) = 1.74123 ± 0.17982 B 

The shape of the Electroweak fit A% 2 around its global minimum follows a Gaussian distri- 
bution. Therefore, it is modelled by a single observable log 10 (m#-o/GeV), linearly correlated to 
mf ie and Aaf^ dl while the two remaining parameters, mz and a s [rriz\ are considered constant. 



The linear correlation factors r[log 10 (m#o), Aa 



(5) 
hadi 



0.261 and r[log 10 (m#o), m 



polei 



-0.494 

are derived from the shape of the multidimensional Ax 2 in the neighbourhood of its global 
minimum, as: 



C[x,y) 



dy* a x 

dx <7„ 



(51) 



where y*[x] is the parameter value that minimises the Ax 2 for the fixed value of x = log 10 (m#o/GeV) 
and where the derivative is taken at the best guess value x*, the one yielding the global mini- 
mum. 

The left and right couplings, g\ and g b R , are expressed as a function of sin 2 6 b eE and the 
Veltman's parameter pb : 



9 b n 



-y/PbQb sin 2 



effy 



(52) 
(53) 



whith T| = 1/2 and Qb = —1/3 the third component of the weak isospin and the electric charge 
of the b quark. The SM parameters Kb, sin 2 (8 b s ) and pb are thus parameterized by polynomials 



and Aa[2j! in order 



of the 3 relevant input parameters of the Electroweak fit mf° e , log 10 (mfl-o 
to reproduce ZFITTER results [HI]. Variations with rnP° le and Aa[ 5 a ' d are linear up to 5a from 
their expectation. Variations with log 10 (m#o) could however only be reproduced by a 4 th order 
polynomial. The parameterization is written as: 




.4 



Z^bb 



f A t \ 
A 
A 



A 2 



1.333120 N 
| 0.232641 
.0.994306, 



(54) 



where Ax are the variations of the parameters m^ ole , Aa^ d and log 10 (m#o) (denoted X = 
t,a,H), around their central values (respectively 172.4 GeV, 0.02758 and 1.74123). A z ^ b l ls 
a 3 x 6 matrix whose coefficients are given in the transposed matrix of eq. [551 With these 
numerical values one reproduces the predictions for K b = 1.33312, sin 2 b cS = 0.232641 and 
Pb = 0.994306 to a relative accuracy of 10~ 6 , for values of the input parameters within a 3cr 
range from their central values. 



4 It has been checked that the global EW fit results 
into account. 



were accurately reproduced when taking all inputs 
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x 1 


x 2 




3 

X 




4 

X 




Pi 


0.953342 


-0.197937 


5.50096 • 


lO" 2 


-3.22697 


10~ 3 


1.30788- 


10" 4 


Q2 


-0.219441 


-2.72075 ■ 10" 2 


-5.26861 


• lO" 3 


-2.74128 


10" 4 


-7.53345 


• 10" 6 



Table 5: Polynomial coefficients for the parameterization of the 2HDM(II) prediction of R^. Coeffi- 
cients are listed in increasing power of x. 



A 1 



( 


1.97262 ■ 


lO" 4 


-6.07684 • 10" 6 


-7.34121 


io- 5 \ 




-0.6* 


35225 


0.348494 


1.73585 


lO" 3 




-1.32964- 


lO" 3 


8.52074 • lO" 4 


6.36640 


lO" 4 




-3.91140 ■ 


lO" 3 


1.57997- 10~ 4 


-2.59313 


lO" 3 




9.57438 ■ 


10- 4 





3.99581 


10- 4 


\ 


2.62977- 


10" 4 





2.56867 


10-7 



(55) 



C.2 Charged Higgs contribution 

The charged Higgs contributions can be implemented through a redefinition of the coupling 
constants, as given in Eq. (1271) . where the function F z , at 2 loops is: 



F z [x] = m + ^^f2[x}. 



(56) 



with x = m 2 (fi )/m 2 H+ . Computations are performed at the scale /i = 2m w as for b — > sj, 
in Appendix B. The two loops correction, f 2 , was computed in ref. [5U]. To avoid di-logarithm 
computations, we parameterized the corresponding terms in f 2 by rational fractions, as: 

x 2 — x — x ln(x) 



/iN 



[X 



\2fx[x\ -6x ( m t (fi ) 2 
/aN = : In 



x 



/'o 



(57) 
(58) 



where P 2 and Q 2 are 4 th order polynomials whose coefficients are given in Table [5] The latter 
parameterization reproduces the results of [50] with an accuracy better than 0.5% for x in a 
range from 10 -5 to 10 5 . 
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